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Kinetic theory of two-dimensional point vortices with collective effects 
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We develop a kinetic theory of point vortices in two-dimensional hydrodynamics taking collective 
effects into account. We first recall the approach of Dubin & O'Neil [Phys. Rev. Lett. 60, 1286 
(1988)] that leads to a Lenard-Balescu-type kinetic equation for axisymmetric flows. When collective 
effects are neglected, it reduces to the Landau-type kinetic equation obtained independently in our 
previous papers [PH. Chavanis, Phys. Rev. E 64, 026309 (2001); Physica A 387, 1123 (2008)]. We 
also consider the relaxation of a test vortex in a "sea" (bath) of field vortices. Its stochastic motion 
is described in terms of a Fokker-Planck equation. We determine the diffusion coefficient and the 
drift term by explicitly calculating the first and second order moments of the radial displacement of 
the test vortex from its equations of motion, taking collective effects into account. This generalizes 
the expressions obtained in our previous works. We discuss the scaling with A'^ of the relaxation 
time for the system as a whole and for a test vortex in a bath. 



r-| ■ I. INTRODUCTION 

o ; 

^ ' The point vortex model |l| is interesting not only because it provides a simplified model of two-dimensional (2D) 

hydrodynamics with potential applications to large-scale geophysical and astrophysical flows, but also because it 

constitutes a fundamental example of systems with long-range interactions. Systems with long-range interactions 

are numerous in nature (self-gravitating systems, non-neutral plasmas, two-dimensional vortices,...) and they present 

C/3 . striking features that are absent in systems with short-range interactions [2t3| • Their dynamics and thermodynamics 

are actively studied at present Q . The statistical mechanics of the point vortex ga s wa s first considered by Onsager 

Q who discovered that negative temperature states are possible for this system J102J ]. At negative temperatures 

(corresponding to high energies) , point vortices of the same sign have the tendency to cluster into "supervortices" , 

n^ ] similar to the large-scale vortices (e.g. Jupiter's great red spot) observed in the atmosphere of giant planets. The 

r^ ■ qualitative arguments of Onsager [6[ were developed more quantitatively in a mean field approximation by Joyce & 

Q \ Montgomery [^, Q , Kida [1^ and Pointin & Lundgren [ll|, [l2] ; and by Onsager himself in unpublished notes [13 ■ 

O . The statistical theory predicts that the point vortex gas should relax towards an equilibrium state described by the 

Boltzmann distribution. Specifically, the equilibrium stream function is solution of a Boltzmann-Poisson equation. 

,_^ , Many mathematical works |14l4l8l ] have shown how a proper thermodynamic limit can be rigorously defined for the 

^ ■ point vortex gas (in the Onsager picture). It is shown that the mean field approximation becomes exact in the limit 

C \ N -^ +00 with 7 ^ 1/N and _R ~' 1 (where N is the number of point vortices, 7 is the circulation of a point vortex and 

\l ' R is the system size). However, the equilibrium statistical theory is based on some assumptions of ergodicity and it is 

'^ \ not proven that the point vortex gas is ergodic [l^- Therefore, it is not clearly established that it will relax towards 

. . the Boltzmann distribution for t — > +cxj. To settle this issue, and to determine the relaxation time (in particular its 

C^^ ' scaling with N), we must develop a kinetic theory of point vortices. 

^^ . A kinetic theory has been developed by Dubin & O'Neil [20| in the context of non-neutral plasmas under a strong 

, ' magnetic field, a system isomorphic to the point vortex gas. They consider axisymmetric mean flows and take 

. . , collective effects into account. Their approach is formal (but powerful) and uses the methods developed in plasma 

^ ■ physics. They perform a quasilinear theory of the Klimontovich equation, introduce Fourier-Laplace transforms and 

work in the complex plane. They obtain a kinetic equation (|^T|) that is the counterpart of the Lcnard-Balescu [2l|, ^^ 

equation in plasma physics. Their work has been continued in [23|-l26|, albeit in a different perspective. 

C^ ' Independently, we have developed a kinetic theory of point vortices |27h31| by using an analogy with the kinetic 

theory of stellar systems |32| - [3a l . The originality of our approach is to remain in physical space so as to treat flows 

that are not necessarily axisymmetric. We obtained a general kinetic equation (|44p that applies to arbitrary mean 

flows. We derived it by different methods such as the projection operator formalism [25, the BBGKY hierarchy [2^ 

and the quasilinear theory [23 . Interestingly, the structure of this kinetic equation bears a clear physical meaning in 

terms of generalized Kubo relations. This equation is valid at the order 1/A^ and, for A^ -^ -|-cxd, it reduces to the 

(smooth) 2D Euler equation (fT5|) which describes the collisionless evolution of the point vortex gas. For axisymmetric 

mean flows, we obtained a simpler and more explicit kinetic equation (|42[) that is the counterpart of the Landau |37| 

equation in plasma physics. Our formalism |29l . l3l| is simple and clearly shows the structure of the kinetic equation 

and the physical origin of its different components. Furthermore, it avoids the use of Fourier-Laplace transforms and 

provides a direct and transparent derivation of the kinetic equation without going into the complex plane (which, 

somehow, hides the basic physics). The derived kinetic equation conserves the circulation, the energy and the angular 



momentum. It also increases monotonically the Boltzmann entropy (iJ-theorem) . The colhsional evolution is due to a 
condition of resonance between distant vortices and the relaxation stops when the profile of angular velocity becomes 
monotonic, even if the system has not reached the Boltzmann distribution of statistical equilibrium [33 ■ Since the 
kinetic equation is valid at the order 1/A^, this "kinetic blocking" implies that the relaxation time for an axisymmetric 
distribution of point vortices is larger than -/Vi_D, where t]j is the dynamical time. In fact, there is no proof so far 
from the kinetic theory that the system relaxes towards the Boltzmann statistical equilibrium for t — > +cxd. To settle 
this issue, we need to develop the kinetic theory at higher orders (taking into account more complex three-body, 
four-body... correlations between vortices). By contrast, for a non-axisymmetric evolution, there arc potentially more 
resonances and the relaxation time can be shorter, of the order of the natural timescale Nto [3l|. We also considered 
the stochastic motion of a test vortex in a "sea" (bath) of field vortices |27H3]| . It can be described in terms of a 
Fokker-Planck equation involving a diffusion term and a drift term. For a thermal bath, they are connected to each 
other by an Einstein relation. The diffusion coefficient (resp. drift term) is proportional to the vorticity distribution 
(resp. to the gradient of the vorticity distribution) of the field vortices and inversely proportional to the local shear, 
a feature first noted in [231 ■ The distribution of the test vortex relaxes towards the distribution of the field vortices 
(bath) on a timescale {N /In N)tD- 

A formal limitation of our approach is that it neglects collective effects. This is not so crucial as long as we are just 
interested in the main properties of the kinetic equation that do not sensibly depend on collective effects (as we shall 
see, collective effects essentially amount to replacing the "bare" potential of interaction by a "dressed" potential of 
interaction, without altering the overall structure of the kinetic equation). However, collective effects may change the 
results quantitatively |103| . In this paper, we generalize our kinetic theory so as to take collective effects into account. 
We also extend the kinetic theory to a general class of potentials of interaction between point vortices (in particular, 
they may include screening effects that are important in geophysics [40|). In Sec. IIIIl we consider the evolution of 
the system as a whole and re-derive the Lcnard-Balcscu-type kinetic equation obtained in Ref. [20] . This is done for 
the sake of completeness and in order to obtain intermediate results that will be useful in the sequel. We also show 
the connection with the Landau-type kinetic equation that we independently derived in Refs. [27H31J with another 
method when collective effects are neglected. In Sec. IIVI we consider the relaxation of a test vortex in a bath of field 
vortices and derive a Fokker-Planck equation (this Fokker-Planck approach is not considered in [20]). We determine 
the diffusion coefficient and the drift term by explicitly calculating the first and second order moments of the radial 
displacement of the test vortex from its equations of motion, taking collective effects into account. This generalizes 
the expressions obtained in our previous works [27K31[ where collective effects are neglected. Finally, we show the 
connection between the Fokker-Planck equation describing the relaxation of a test vortex in a bath and the kinetic 
equation describing the evolution of the system as a whole and compare the relaxation time in these two situations. 

There exists many analogies between systems with long-range interactions although they are o f a v ery different 



nature. We have already mentioned analogies between plasmas, stellar systems and point vortices [10J|. There also 
exists analogies with a toy model of systems with long-range interactions called the Hamiltonian Mean Field (HMF) 
model [4j] . The kinetic theory of this model is developed in [3, [3l|, l44| - l5g | and it presents features that are similar to 
those observed earlier for point vortices [20, [23, [23] ■ In particular, the Lenard-Balescu collision term vanishes in one 
dimension due to the absence of resonances (a result known long ago in plasma physics [5ll |52| and rediscovered in 
the context of the HMF model [4^, [43] ) so that the relaxation time of the spatially homogeneous HMF model is larger 
than Nto (it can be reduced to Nto if the system is spatially inhomogeneous due to additional resonances [3ll. [53j). 
On the other hand, the relaxation of a test particle in a bath is described by a Fokker-Planck equation involving a 
term of diffusion and a term of friction. This is similar to the Fokker-Planck equation derived by Chandrasekhar 
[33l . [34| in stellar dynamics and by Chavanis |27[ - [31| for the point vortex system. This is also a particular case (in one 
dimension and for a potential of interaction truncated to one Fourier mode) of the general Fokker-Planck equation 
derived by Hubbard [54| in plasma physics. In order to stress the analogies between these different topics, we shall 
use a presentation that is similar to the one given in the review of Campa et al. [5[ on long-range interacting systems. 

II. THE POINT VORTEX MODEL 

Let us consider a system of point vortices with individual circulation 7 moving on the infinite plane. Their dynamics 
is described by the Hamilton-Kirchhoff Jj, _55i] equations 

dxi dH dyi dH 

dt dyi ' dt dxi ' 

with the Hamiltonian 

H = ^^Y.<\v,^v,\). (2) 

i<j 



Wc note that the coordinates x and y of the point vortices are canonically conjugate. The standard potential of 
interaction between point vortices is solution of the Poisson equation Au = — (S. In an infinite domain, it is given 
by u(|r — r'l) = — (l/27r)ln|r — r'|. However, we shall let the function u arbitrary in order to treat more general 
situations. For example, in the quasigeostrophic (QG) approximation of geophysical fluid dynamics, the Poisson 
equation is replaced by the screened Poisson equation Au — k\u = —6 where fc^ is the Rossby radius [40|. Equation 
([1]) can be rewritten in vectorial form 

dv ■ 
7^ = -z X V.i/, (3) 

where z is a unit vector normal to the plane on which the vortices move. In an infinite domain, the point vortex 
system conserves the total circulation P = N^, the energy E ~ H and the angular momentum L = ^^ jrf (the center 
of "mass" R = X^i T'"* ^^ ^^^'^ conserved and taken as the origin O of the system of coordinates). We introduce the 
discrete vorticity field 



^d 



(r,i)=7E'5('--^^W)- (4) 



Differentiating this expression with respect to time and using the equation of motion ([3]) , we get 

diOd 



dt 



^(zx V,i/)-V(5(r-r,(t)). (5) 

i 

This equation can be rewritten 

dtod 
dt 



^ = V • ^(z X W.,H)6{r - r.,{t)) = 7V • ^^(z x V^,)<5(r - r,(i)), (6) 



where 

Mr,t)=Jui\r~r'\)iOd{r\t)dr\ (7) 

is the discrete stream function induced by the discrete vorticity field (we have used the J-function to replace 
VV'd(ri(i),i) by V?/'d(r,i) in the last equality of Eq. ©). Introducing the discrete velocity field Ud = — z x Vtpd, we 
obtain 

^ + ViuJdUd) = 0. (8) 

ot 

Finally, using the incompressibility condition V • u^; = 0, we can rewrite this equation in the form 

^ + Ud • Vcjd = 0. (9) 

This equation is exact and contains the same information as the Hamiltonian system ((T])-([2]). This is the counterpart 
of the Klimontovich equation in plasma physics [56[ . This is also the 2D Euler equation from which the point vortex 
model is issued [1| . 

We now introduce a smooth vorticity field a;(r, i) = ((jjd(r, i)) corresponding to an average of ajd(r,t) over a large 
number of initial conditions. We then write w^ = w + ^w where (5w denotes the fluctuations around the smooth 
distribution. Similarly, we write ipd = i-' + Sip and u^ = u + Su. Substituting these decompositions in Eq. (|9]), we 
obtain 

duj dSuj ^ „^ ^ „ ^ „^ „ ,,„, 

— - + — h u • Vw + u • VSuj + Su-Vuj + Su- VSuj = 0. (10) 

ot dt 

Taking the average of this equation over the initial conditions, we get 

-^+u-\/uj ^ -V ■ (SljSu), (11) 

ot 

where the right hand side can be interpreted as a "collision" term. Of course, it does not correspond to direct 
collisions between vortices but to distant collisions, or correlations, due to the long-range interaction. Subtracting 
this expression from Eq. (|10|). we obtain 

or 

-^ + u • V(5w + (5u • Vw = V • (SujSu) - V • {6uj6u). (12) 



These equations are still exact. We now consider the thermodynamic limit A^ -^ +00 with 7 ^ 1/N. In this limit, we 
have Soj ^ 1/VN and i5'0 ~ ^/VN and we see that the collision term scales like 1/A^. Therefore, for N -^ +00, Eq. 
PT|) reduces to the (smooth) 2D Euler equation 

^ + u-Vw = 0, u = -zxVip, ip{r,t)= u{\r-r'\)uj{r',t)dr'. (13) 

The 2D Euler equation is obtained when "collisions" (more properly granular effects, discreteness effects or finite N 
corrections) between point vortices are neglected. The 2D Euler equation can also be obtained from the BBGKY 
hierarchy by neglecting correlations between point vortices [29|. This is valid, for a fixed interval of time, when 
N — > +00. In that case, the mean field approximation becomes exact and the iV-body distribution function can 
be factorized as a product of N one-body distributions, resulting in Eq. ([T^. The 2D Euler equation ([T^ is the 
counterpart of the Vlasov equation in plasma physics j56j and stellar dynamics [57| . Starting from an unsteady or 
unstable initial condition, the 2D Euler equation is known to develop a complicated mixing process leading to the 
formation of a quasi stationary state (QSS) on a coarse-grained scale |58l . |59[. This takes place on a very short 
time scale, of the order of the dynamical time tjj ^ (^~^, where w ^ T / E? is the typical vorticity (i? is the system size 
[l05l |). This process, which is purely coUisionless and driven by mean field effects, is known as violent relaxation. A 
statistical theory of the 2D Euler equation has been developed by Miller [60| and Robert & Sommeria [6l| to predict 
the QSS that results from violent relaxation assuming ergodicity. This is the counterpart of the statistical theory of 
the Vlasov equation introduced by Lynden-Bell [G^l in astrophysics to describe the violent relaxation of coUisionless 
stellar systems. The process of violent relaxation towards a non-Boltzmannian QSS has also been exte nsiv ely studied 
for the HMF model p. In some cases, violent relaxation is incomplete due to lack of efficient mixing [106| . 

On longer timescales, "collisions" (more properly correlations) between point vortices develop and the system 
deviates from the 2D Euler dynamics. For t — >■ -|-c», we expect that the system will reach a statistical e quil ibrium 
state. In the thermodynamic limit N — > +00 with 7 ^ 1/-/V, the mean field approximation is exact [14l4l8j . The 
statistical equilibrium state in the microcanonical ensemble is obtained by maximizing the Boltzmann entropy S = 
(w/7) ln(w/7) dr at fixed circulation F = J w dr, energy E = (1/2) J ujij} dr and angular momentum L = J ujr'^ dr 

|38|. This variational principle determines the most probable macrostate w(r). Writing the first order variations as 
6S — f36E — a6r— {\/2)PQ,l5L = where /?, a and $7^ are Lagrange multipliers, we obtain the Boltzmann distribution 

Lc = Ae~f^^-^', i;{r)= [ u{\r-r'\)ujir')dr', (14) 



where V'* = V' + (1/2)^l''^ is the relative stream function taking into account the invariancc by rotation of the system. 
For the ordinary inte raction between point vortices, the equilibrium stream function is solution of the Boltzmann- 
Poisson equation [107l | . We must remember, however, that the statistical theory is based on an assumption of ergodicity 
and on the postulate that all the accessible microstates are equiprobable. There is no guarantee that this is true for 
the point vortex model [1^ . If we want to prove that the point vortex gas truly reaches the Boltzmann distribution 
(|14p. and if we want to determine the relaxation time (in particular its scaling with the number A^ of point vortices), 
we must develop a kinetic theory of point vortices. 

III. COLLISIONAL EVOLUTION OF THE SYSTEM AS A WHOLE 

A. Lenard-Balescu and Landau-type kinetic equations 

As we have previously indicated, the "collisions" between point vortices can be neglected for times shorter than 
Nto where t^ is the dynamical time. If we want to describe the evolution of the system on a longer time scale, we 
must take finite N corrections into account. At the level 1/N, we can neglect the quadratic term on the right hand 
side of Eq. (|12p . Indeed, the left hand side is of order 1/\/N and the right hand side is of order 1/A^ <C 1/Vn. We 
therefore obtain the set of coupled equations 

^ -h u • Vcj = - V • (SljSu) , (15) 

at 

— ^ + u-V(5u; + (5u-Vw = 0. (16) 

ot 



They form the starting point of the quasihnear theory. If we restrict ourselves to axisymmctric mean flows, we can 
write u(r, i) = u{r,t)eg with u{r,t) = —dip/dr{r,t) = il.{r,t)r, where ^{r,t) = r~^ L u!{r' ,t)r' dr' is the angular 
velocity field (see Sec. 6.1 of JSM])- In that case, Eqs. (fT5|) -(|16 |) reduce to the coupled equations 

duj Id/ d5^\ 

We shall assume that the fluctuations evolve rapidly compared to the transport time scale, so that time variation of w, 
■0 and J7 can be neglected in the calculation of the collision term (Bogoliubov ansatz). For brevity, we shall omit the 
variable t in these quantities. We also assume that the axisymmctric distribution always remain Euler stable, so that 
it evolves only under the effect of collisions and not because of dynamical instabilities. In that case, Eqs. (|17p - p^ 
can be solved with the aid of Fourier-Laplace transforms and the collision term can be explicitly calculated. This is 
the kinetic approach developed by Dubin & O'Neil ^^. We shall recall the main steps of their analysis by using a 
close parallel with the derivation of the Lcnard-Balescu equation in plasma physics [5i, ,56,, 67, ,68i] . This will allow us 
to set the notations and derive useful intermediate results that will be needed in the sequel. 
The Fourier-Laplace transform of the fluctuations of the vorticity fleld 6uj is defined by 

dQin,r,a)= — df e-*("''-"*)(5w(^, r, i). (19) 



2^ JO 



This expression for the Laplace transform is valid for Im(cr) sufficiently large. For the remaining part of the complex 
a plane, it is defined by an analytic continuation. The inverse transform is 



v^ / da 



Suj{9,r,t)= }_^ j —e^(-'-^')5G:{n,r,a), (20) 

n— — oo ' 



^ /c2^ 



where the Laplace contour C in the complex a plane must pass above all poles of the integrand. Similar expressions 
hold for the fluctuations of the stream function (5-0. If we take the Fourier-Laplace transform of Eq. (fTS]). we find that 



1 /-I 

— (5w(n, r, 0) — iaSdiin, r, a) + inflSioin, r, a) + in——-6ip{n, r, a) = 0, (21) 

r or 

where the first term is the spatial Fourier transform of the initial value 

Su;{n,r,0)^J — e~^"'Sio{e,r,0). (22) 

The foregoing equation can be rewritten 

duj{n, r, a) = -^-^^-Sipin, r, a) + — — -, (23) 

ml ^ a t[nil — a) 

where the first term on the right hand side corresponds to "collective effects" and the second term is related to the 
initial condition. The fluctuations of the stream function are related to the fluctuations of the vorticity by an equation 
of the form LSip = —Suj equivalent to Sip = u * Suj (see Appendix E^ . Taking the Fourier-Laplace transform of this 
equation, we obtain C6ip{n,r,a) — -~5u]{n,r,a). Using Eq. (1231) . it can be rewritten in the form |108l | 



C- 



n-TT 

r or 

nQ — a 



rT/ \ (5w(n, r, 0) 

(5V(n,^fT) = -— ^- '-. (24) 

i[n\l — a) 



Therefore, the Fourier-Laplace transform of the fluctuations of the stream function is related to the initial condition 
by 

Hin.r.a)^ T" 2^r' dr' G(n, r, r', a)^^%^, (25) 

Jo i[n\V - a) 



where the Green function (^(n, r, r' , a) is defined by 



C- 



r or 

iQ — a 



G{n,r,r',(T) = - 



S{r - r') 
27rr 



(26) 



The Fourier-Laplace transform of the fluctuations of the vorticity is then given by Eq. (P^ with Eq. ([^5)) . In these 
equations, we have noted J7 for n(r) and Q' for 0(7"'). Similarly, we shall note u! for Lj{r) and cj' for w(r'). To avoid 
confusion, the derivatives with respect to r will be noted d /dr. 

We can use these expressions to compute the coUisional term appearing on the right hand side of Eq. (fTT]). One 
has 



, d5-<p\ 



EE 



da f da' 



2tt r 27r 



■^,^^ine-at)^^{n'e-a't) ^^^(^^ ^^ a)54>{n',r, a')). 



(27) 



Using Eq. dSS]), we find that 



1 9w 



{Suj(n,r,a)S^p(n' ,r,a')) ~ J' {Sip(n,r,a)Sip{n ,r,a')) 

ml — a 



{6u;{n,r,0)6i>{n',r,a')) 
i(nfl — a) 



(28) 



The first term corresponds to the self-correlation of the stream function, while the second term corresponds to the 
correlations between the fluctuations of the stream function and of the vorticity at time t = 0. Let us consider these 
two terms separately. 

From Eq. (|25p . we obtain 



{Sip{n, r, a)6'ip{n' , r, a')) 



+ 00 



„ ,JSLj{n,r',0)Sw{n',r",0)) 



27Tr'dr' 27rr"dr"G{n,r,r',a)G{n',r,r",a')- „ 

Jo [nil' - a)[n'Q" - a') 



Using the expression of the auto-correlation of the fluctuations at t = (see Appendix |B]) given by 

{du;{n, r', 0)<5c2;(n', r" , 0)) = jS^^_JX^-^ujir'), 



we find that 

{Sip{n,r,a)S^{n' ,r,a')) :^ 2tt^ 6n.-r 



27rr' 



+ CXD / /\ 

r' dr' G{n, r, r' , a)G{-n, r, r', cr')-^-— 



{nfl' - a){nn' + a')' 



(29) 



(30) 



(31) 



Considering only the contributions that do not decay in time, it can be shown [Sl, [5a] that [{nQ' — a){nri' + cr')]~^ 
can be substituted by (2'K)^5{a + a')5{a — nfl'). Then, using the property G{—n, r, r' , —a) ~ G{n, r, r' , a)* [2^, one 
finds that the correlations of the fluctuations of the stream function are given by 



((5V'(n,r,cr)JV(n',r,cr')) = (271)^7 (5„._„/(5(cr + cr') / r' dr' \G{n,r,r' ,a)Y5{a - nVt')uj{r') 



(32) 



This is an important relation that we wish to emphasize. Similarly, one finds that the second term on the right hand 
side of Eq. (PSJ) is given by 



{6Lo{n,r,Q)5i){n',r,a')) 



^ (27r)^7 5n~n'5{a + a')G{-n, r, r, -a)5{a - nn)uj{r). 



i(nn — a) 
From Eq. p2p , we get the contribution to ([27| of the first term of Eq. psp . It is given by 

n+oo „2 1 Sti 



Slj 



de / J ^ ' ' ^^Jc27r Jq nQ-a 



\G{n,r,r',a)\H{a - nn')uj{r'). 



Using the Plemelj formula, 



lim 

e->o X — a± iO+ 



V- 



=p i-K5{x — a). 



(33) 



(34) 



(35) 



where V denotes the principal value, integrating over a and using the identity 



5{\x) = ^Mx) 



(36) 



we obtain 



r dr 



Suj^\ ^-2Tr'\y2f r'dr'\n\\G{n,r,r',nn)fS{n-n')uj{r') 
o(^ / I „ io 

From Eq. (|33|) , we get the contribution to ((27)) of the second term of Eq. ([28|) . It is given by 



5lj 



-Q^j ='^'^1^1 ■:^nlTaG{n,r,r,a)5{a -nVL)uj{r). 



(37) 



(38) 



To determine ImG(n,r, r, cr), we multiply Eq. p6p by 2T:r G{—n,r,r' , —a), integrate over r from to +oo and take 
the imaginary part of the resulting expression using the Plemelj formula. This yields 



ImG(n,r,r,a)=T2TT^ / dr' n\G(n,r,r' ,a)\^S(a - nn')—-. 
Jo or 

Substituting this expression in Eq. ([55)1 and using Eq. ([M|) . we obtain 

(^Sco^'j =2n'-jY.£^r'dr'\n\\Gin,r,r',nnrS{n-n')u;^^. 
Finally, regrouping Eqs. ()17p . (|37| and (|40|) . we end up on the kinetic equation 



'dt 



r dr 



2-V^E 



r'dr' \n\\G{n,ry ,nn)\H{n - n') [ -u'^ - ^u^ 

r or r or 



(39) 



(40) 



(41) 



This equation, taking collective effects into account, was derived by Dubin & O'Neil [20]. It is the counterpart of the 
Lenard-Balescu equation in plasma physics. When collective effects are neglected, we independently obtained in Refs. 
[28l - l3ll | a kinetic equation of the form 



'dt 



r dr 



2-V^E 



r' dr' \n\ur,{r, r'fb{n - Q') -uj'- ju—j 

\r or r or 



(42) 



For the usual potential of interaction, the sum over n can be done explicitly (see Appendix \Q and we get 



— - ~ — — — -^ / r' dr' In 
at 4 r or Jq 






5(n - n') i -u}'— - —oj—] 

^ r dr r' dr' ' 



(43) 



where r< (resp. r>) is the min (resp. max) of r and r' . This kinetic equation was derived with a method that does 
not use Fourier-Laplace transforms. It arises from the generalized kinetic equation 



duj N -I 
'dt^ N 



u-\/uj^ —^ I dr I driVil -^ 0)g{t, t - r) 



^(1^0)^— + V"(0^1)— 7 



uj{r,t-T)-{ri,t~T), 

7 



(44) 



derived in [28|, |29| , which is valid for flows that are not necessarily axisymmetric and not necessarily Markovian [109l | . 
Equation (J43|) is the counterpart of the Landau equation in plasma physics. The connection between Eqs. (|4ip and 
dH]) is clear. If we neglect collective effects, Eq. (^5)) reduces to CGbare{n,r,r') = —6{r — r')/2Trr (it does not depend 
on a anymore) so that Gbarein,r,r') is just the Fourier transform of the "bare" potential of interaction u (i.e. the 
Green function of L). We thus have Gbarein.r.r') = Un.{r,r') (see Appendix \K\ . When collective effects are taken 
into account, the "bare" potential Gbare{n,r,r') = Un{r,r') is replaced by a "dressed" potential G{n,r,r',a) without 
changing the overall structure of the kinetic equation. This is similar to the case of plasma physics where the bare 



potential u{k)'^ in the Landau equation is replaced by the dressed potential u(A:)^/|e(k, k- v)|^, including the dielectric 
function, in the Lenard-Balescu equation. In plasma physics, collective effects are important because they account 
for screening effects and regularize, at the scale of the Debye length, the logarithmic divergence that occurs in the 
Landau equation. In our context, there is no divergence in Eq. (j43|) so that collective effects may be less crucial than 
in plasma physics. So far, their influence is not clearly understood. 
Finally, we note that the kinetic equation can be rewritten 

^ = 27rV|-y^ r'dr'x(r,r',i)<5Mr,i)-f](r',i))(^-|---^jc.(r,iV(r',i), (45) 



with the notation 



X{r,r',t) = Y, \n\\Gin,r,r',nn{r,t))\^ (46) 



When collective effects are neglected %('',''', is replaced by Xbareifjr') = J^n \M''^n{r,r')'^. For the usual potential 
of interaction (see Appendix E| . it reduces to Xbareir, r') = — (I/Stt^) In [l — (r</r>)^] . 

B. The relaxation time of the system as a whole 

The kinetic equation (|^5|) is valid at the order 1/N so that it describes the "coUisional" evolution of the system on 
a timescale ~ Ntu. This kinetic equation conserves the total circulation F, the energy E and the angular momentum 
L. It also monotonically increases the Boltzmann entropy S (H-thcoveia). The Boltzmann distribution (|14p is always 
a steady state of this kinetic equation, but it is not the only one: any vorticity distribution which is associated 
with a monotonic profile of angular velocity is a steady state of Eq. (|45p . The kinetic equation admits therefore an 
infinite number of steady states! As explained in previous works, the coUisional evolution of the point vortex gas 
(at the order 1/A^) is due to a condition of resonance between distant orbits of the point vortices. For axisymmetric 
systems, the condition of resonance, encapsulated in the i5-function, corresponds to i^{r\t) = 0(r, i) with r' ^r (the 
self-interaction at r' = r does not produce transport since the term in parenthesis in Eq. (|45p vanishes identically). 
The vorticity profile u>{r,t) at position r will change under the effects of "collisions" if there exists point vortices 
at r' 7^ r which rotate with the same angu lar v elocity as point vortices located at r. This is possible only if the 
profile of angular velocity is non- monotonic [llOj. Therefore, the evolution is truly due to distant collisions between 
point vortices (this is strikingly different from the case of plasmas and stellar systems where the collisions are more 
"local"). The evolution stops when the profile of angular velocity becomes monotonic (so that there is no resonance) 
even if the system has not reached the statistical equilibrium state given by the Boltzmann distribution. This "kinetic 
blocking" has been illustrated numerically in [3^. Distant collisions between point vortices have the tendency to 
create a monotonic profile of angular velocity. If the initial condition already has a monotonic profile of angular 
velocity, the collision term vanishes (at the order 1/N) because there is no resonance. The kinetic equation reduces 
to 

so that the vorticity does not evolve at all on a timescale ^ Ntjy. This implies that, for an axisymmetric distribution 
of point vortices, the relaxation time is larger than Ntij. We therefore expect that 

tji > Nt]j (axisymmetric flows) (48) 

Since the relaxation process is due to more complex correlations between point vortices, we have to develop the kinetic 
theory at higher orders (taking into account three-body, four-body,... correlation functions) in order to obtain the 
relaxation time. If the collision term does not vanish at the next order of the expansion, the kinetic theory would 
imply a relaxation time of the order N^tu |3l| . However, the problem could be more complicated and yield a larger 
relaxation time like e^to- In fact, it is not even granted that the system will ever relax towards statistical equilibrium 
(the point vortex gas may be non-ergodic [19{). By contrast, for non-axisymmetric flows, since there are potentially 
more resonances between point vortices [see the complicated kinetic equation (HU], the relaxation time could be 
reduced and approach the natural scaling 

tfi ^ Ntn (non — axisymmetric flows) (49) 

predicted by the first order kinetic theory. This linear scaling has been observed for 2D point vortices with non- 



axisymmetric distribution in |71| . However, very little is known concerning the properties of Eq. (|44p and its 



convergence (or not) towards the Boltzmann distribution. It could approach the Boltzmann distribution (since 
entropy increases), without reaching it exactly. 

Similar results have been found for one-dimensional plasmas and stellar systems and for the HMF model (see 
the conclusion for more details). For spatially homogeneous one dimensional systems, the Lenard-Balescu collision 
term vanishes (no resonance) so that the relaxation time is larger than Nto- For spatially inhomogeneous one 
dimensional systems, since there are potentially more resonances (this can be seen by using angle-action variables 
[53|), the relaxation time can be reduced and approach the natural scaling Nto- For d-dimensional plasmas and 
stellar systems, with d > 1, there arc always resonances since the condition is k • v = k • v'. The Lenard-Balescu 
collision term vanishes only for the Boltzmann distribution and the relaxation time is tfl ^ Ntn (or 1^ ^ {N/ In N)tD 
for stellar systems) corresponding to the first order of the kinetic theory. 

IV. STOCHASTIC PROCESS OF A TEST VORTEX: DIFFUSION AND DRIFT 

A. The Fokker-Planck equation 

In the previous section, we have studied the evolution of a system of point vortices as a whole. We now consider 
the relaxation of a test vortex is a "sea" of field vortices with a steady axisymmetric vorticity profile uj{r). The 
field vortices play the role of a bath. We assume that the field vortices are either at statistical equilibrium with 
the Boltzmann distribution (thermal bath), in which case their distribution does not change at all, or in a stable 
axisymmetric steady state of the 2D Euler equation with a monotonic profile of angular velocity (as wc have just seen, 
this profile does not change on a timescale of order Ntn)- We assume that the test vortex is initially located at a radial 
distance rg and we study how it progressively acquires the distribution of the field vortices due to distant "collisions" . 
As we shall see, the test vortex has a stochastic motion and the evolution of the distribution function P{r,t), the 
probability density of finding the test vortex at radial position r at time f , is governed by a Fokker-Planck equation 
involving a diffusion term and a drift term that can be analytically obtained. The Fokker-Planck equation may then 
be solved with the initial condition P{r^ 0) = d{r — ro)/27rr to yield P{r, t). This problem has been investigated in our 
previous papers |27H3l| , but we shall give here a direct and more rigorous derivation of the coefficients or diffusion 
and drift, taking collective effects into account. 

The equations of motion of the test vortex are 

dr IdSi' d0 ^, , IdS^ 

They include the effect of the mean field fl{r) which produces a zeroth-ordcr net rotation plus a stochastic component 
(5-0 of order I/Vn (fluctuations) which takes into account the deviations from the mean field. They can be formally 
integrated as 

r{t)^r+ r dt' ^^{r{t'),9{t'),t'), (51) 



r{t') de 

e{i) =(^ + [ dt' nir{t')) - 1^ dt' -l^^(r(t'), 0{t'),t'), (52) 

where we have assumed that, initially, the test vortex is at {r,d). Note that the "initial" time considered here does 
not necessarily coincide with the original time mentioned above. Since the fluctuation (5'0 of the stream function is a 
small quantity, the foregoing equations can be solved iteratively. At the order 1/A^, which corresponds to quadratic 
order in Sip, we get for the radial distance 



in iiM\ 



ii f i," L^l^fr.e + at'.t'&..e + at".n + f dt' f cit"\?^(r.e + nt'.t'&:e + Qt",t") 



JO ' '^'^ ^'' Jo Jo ^^ drd9 89 

dt' f dt"\^ir,9 + nt',t')^{r,9 + nt",t"). (53) 



r 



3 d9 ' ' ' ' 89 
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As the changes in the radial distance are small, the dynamics of the test vortex can be represented by a stochastic 
process governed by a Fokker-Planck equation [T^l- If we denote by P{r, t) the density probability of finding the test 
vortex at radial distance r at time f, normalized such that /^ P{r, t)27Tr dr = 1, the general form of this equation is 



dt r dr \ dr 



The diffusion coefficient and the drift term are given by 



1 d 
--^{rPA). 
r or 



(54) 



D{r)= \v^ l{r{t)-rf), 

t^ + oo Zt 



(55) 



Mr) 



lim -(r(t) 



(56) 



In writing these limits, we have implicitly assumed that the time t is long with respect to the fluctuation time but 
short with respect to the relaxation time (of order NtD), so that the expression ((5^ can be used to evaluate Eqs. 
([55)1 and ([55]) . As shown in our previous papers |27H3l| , it is relevant to rewrite the Fokker-Planck equation in the 
alternative form 



dP 

'at 



ld_ 
r dr 



dP 



The total drift term is 



, ^^ dD 

A = V drift H J—, 

dr 



(57) 



(58) 



where Vdrift is the "essential" part of the drift term, while the second term is due to the spatial variations of the 
diffusion coefficient. As we shall see, this decomposition arises naturally in the following analysis. The two expressions 
([S3]) and ([ST]) have their own interest. The expression ([M)) where the diffusion coefficient is placed after the second 
derivative d^{DP) involves the total drift A and the expression (|57l) where the diffusion coefficient is placed between 
the derivatives dDdP isolates the essential part of the drift Vdrift- We shall see in Sec. IIV El that this second form is 
directly related to the kinetic equation (j45|) . It has therefore a clear physical interpretation. 

We shall now calculate the diffusion coefficient and the drift term from Eqs. ([55]) and ([55)) . using the results of Sec. 
nil Al that allow to take collective effects into account. Note that a similar calculation of these terms, directly from the 
equations of motion of a test vortex, has been made in Appendix C of [33|, neglecting collective effects. The diffusion 
coefficient can also be obtained from the Kubo formula and the drift term from a linear response theory |27H29l | . 



B. The diffusion coefficient 



We first compute the diffusion coefficient defined by Eq. ([55|). Using Eq. ([53)) . we see that it is given, at the order 

1/A^, by 



D='rdt'rdt"i'-^i^'^'-^'^'^''^' 



2rH 



\ 36 



By the inverse Fourier-Laplace transform, we have 



.{r,e + m',t')-^{r,e + nt",t 



(59) 



l^{r,e + nt', t')^{r, e + m", t")\ =-EE/^/^ ^^'^^nie+nt')^^^.t'^^n'iO+nt")^^^a't" 

x{STpin,r,a)d^l!{n',r,a')). (60) 

Substituting Eq. (1321) in Eq. (j60p . and carrying out the summation over n' and the integrals over a' and a, we end 
up with the result 

/^^r,0 + nt\t')^{r,e + nt",t")\^2njy r°°r'dr'n2e™(^-"')(*'-*")|G(n,r,/,n^')l'^(^')- (61) 
\ de 80 I ^Jo 



11 



This expression shows that the correlation function appearing in Eq. ((6T|) only depends on \t' — t"\. Using the identity 

t nt nt nt nt 

dt' dt" f{\t' -t"\) = 2 dt' dt" f{\t' ^t"\)^2 ds{t^s)f{s), (62) 



(63) 



/o Jo Jo Jo Jo 

we find, for t — > +oo, that 



-^r(f'^'^'°'t<— ')- 



(64) 



This is the Kubo formula for our problem. Replacing the correlation function by its expression (|6ip. we get 

1 p + ca /• + 00 

D = 2n-f^ dsY r'dr'n^e"'^^^-^^'^''\G(n,r,r',nn')\^u;{r'). 

r^ Jo ^Jo 

Using the symmetry s — )■ — s, we can replace /„ ds by (1/2) J_ ds. Then, using the identity 

/■+°° df 



and Eq. p6p . we obtain the final expression 

/O 



1 f~^°° 

D^2TT^j—y2 r'dr'\n\\G{n,r,r,nn)\^S{n-n')uj{r'). (66) 

''' „ JO 



C. The essential part of the drift term 

We now compute the drift term defined by Eq. ([55]) . We need to keep terms up to order 1/A^. From Eq. ([55]) . the 
first term to compute is 

By the inverse Fourier-Laplace transform, we have 

^ ^ (r, + l^f, t')) = ^ E ^ ^ ,e™(^+"*')e--' (<5^(n, r, a)) . (68) 

Using Eq. (^5]) . we find that 

{di^{n,r,a)) ^ f 2^r' dr' ■^^^%^^^G(n, r, /, ct). (69) 

Jo i(nS2'-cr) 

Now, using the fact that the test vortex is initially located in (r, 6), so that {6uj{9' , r' , 0)) = 7(5(6*' — 6)5{r' — r)/r, we 
obtain from Eq. (f^ . the result 

(to(n,r',0))=7e-™'^^^-^. (70) 

27rr 

Substituting these expressions in Eq. (|68| . we get 

Therefore, the drift term ([ST]) is given by 

^^-^4 f dt' Y. f ^ne^(-^-^^*'G{n,r,r,a)-^. (72) 

'"* Jo „Jc^T^ nn-a 
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Wc now use the Plcmclj formula to evaluate the integral over a. The term corresponding to the imaginary part in 
the Plemelj formula is 

Af^^m-I- dt'Y, — ne*^"""'')*'G(n,r,r,cr)(5(nf]-a). (73) 

Integrating over a and i', we obtain 

A^'') ^_JL^ „IinG(n,r,r,nrj). (74) 

n 

The term eorresponding to the real part m the Plemelj formula is 
Integrating over t', we can convert this expression to the form 

/ + CXj 7 -| -I 
^nG{n,r,r,a)—- {z sin [(nH - c7)t] + eos [(nf^ - a)f] - 1} . 

For t — >■ +00, using the identity (as in Appendix C of [38|): 



and integrating over a, we obtain 



D. The part of the drift term due to the spatial inhomogeneity of the diffusion coefficient 

In the evaluation of the total drift, at the order 1/A^, the seeond term to compute is 



ir-r-r-^f^(^<— ''-'''^^ 



An^- dt' dt" dt'"-—{^^{r,9 + nt',t')^ir,e + nt"\n). (80) 

'0 "'0 "'0 



(76) 



1 — cos(te) -, , ,„„, 



Af'> =-—"^nlmGin,r,r,nn), (78) 



which is the same as Eq. ([74|) . Finally, summing Eqs. fM]) and ([78| and using Eq. ([39| . we find that 

1 f + OO J / 

Ai^2-/TT^-Y dr'\n\\Gin,r,r',nn)fS{n-n')—-. (79) 

r '^ Jo dr' 

As we shall see, it corresponds to the "essential" part of the drift term denoted Vdrift hi Eq. 



^2^1 ^i7i(e+nt') ^-iat' ^in' (9+nt'") ^-la' t'" 



By the inverse Fourier-Laplace transform, we have 

' n n' ^ ^ 

x{Srp{n,r,a)Si}{n',r,a')). (81) 

Substituting Eq. (1321) in Eq. (|8ip, and carrying out the summation over n' and the integrals over a' and ct, we end 
up with the result 



(r,e + nt' ,t')^(r,e + nt'" ,t'")\ ^i2n-iY [ °° r' dr' n^e"''-''-'''^^*'-*"'^\G(n,r,r',nn')\^uj(r'). (82) 
do / V Jo 
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This expression shows that the correlation function appearing in Eq. ((80)) only depends on \t' — 1"'\. Using the identity 

dt" f dt'" f{\t' - r I) = / dt'" it' - t'")f{\t' - t"'|), (83) 

Jo Jo 

we find that 

Au^i2nj- dt' df'Y^ r'dr'{t'-t"')n^^^e'"^''-''^^'~' '>\G{n,r,r',nn')\Mr')- (84) 

i Jo Jo j^ Jo "^ dr 

This can be rewritten 



Aii = 2TT-f- f dt' f df'Y r°° r' dr' n^\^ (e'"^''-'''^^''-'"'A \G{n,r,r' ,nn')\^Lo{r'), 
t Jo Jo „ Jo r^ dr \ J 



(85) 



or, equivalently, 



Aii=2Trj--^ dt' dt"'y^ /dr'n^^e"^"'")^*-* ^\G{n,r,r' ,nQ.')\^uj{r') 

t or Jq Jo „ Jo ''" 

-2n-f-^ C dt' t dt"'J2 r"r'dr'nV"("-"')(*'-*"')|- f^\G{n,r,r',nn')Auj{r'). 



(86) 



Since the integrand only depends on \t' — t"'\, using the identity (|5^. we obtain for t — >■ +oo, 

Ajj =. 2ttj-^ / rfs V / r' dr' n^-^e"'^^-^'>\G{n, r, r' , nrj')l'w(r') 
^^ Jo n Jo ^ 

-27r7 / dsJ2 r'dr'n^e™^^^-^^)^^ I — |G(n,r,r',nrj')nw(r'). (87) 



Using the symmetry s -^> —s, we can replace L ds by (1/2) J_ ds. Then, using the identities (j65p and (j36]), we 
obtain the expression 

A7/==27r27^y / r'dr'\n\—5{n^n')\G{n,r,r',nn')\^Lo{r') 
or ^-^ n r^ 

n ^ 

-2tt''iY1 / r'dr'\n\5{n-n')j- I — |G(n,r,r',nr!')n w(r'). 
In the first term, we recover the diffusion coefficient (j66p so that finally 

Ajj^—- 2tt'jY1 / ^' ^^' l"l'^(" - "')a; [^\G{n,r,r',nn')\'\ iu{r'). 
The third term to compute is 

Ajjj = -^jyt'J^ dt" (^^ir,0 + nt',t')^ir,0 + nt"X)y (90) 



(89) 



This term is just proportional to the diffusion coefficient ([59|) so we get 



2 1 f 

Aiii = D = -4tt^j—Y / r'dr'\n\\G{n,r,r', 

r r* ^^ Jq 



nn)\'^s{n~n')Lo{r'). (91) 



Finally, the fourth and fifth terms to compute are 



Aiv = -^jyt'J^ dt"(^^{r,e + nt',t')^ir,9 + nt",t")), (92) 
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and 

Substituting the inverse Fourier-Laplace transform of the fluctuations of the stream function in these equations, and 
summing the resulting expressions using the fact that the correlation function of the fluctuations of the stream function 
is proportional to J„ _„', we obtain 

Aiv+Av^\- f dt' f df'yy [p /^n2e"(''+"*')e— *'e™'(''+"*")e"-'*"|-(<5V;(n,r,a)5^(n',r,a')). 
r^t Jo Ja „ t^ Jc 27r Jc 27r dr 

(94) 
Substituting Eq. p2p in Eq. (|94p . and carrying out the summation over n' and the integrals over a' and a, we end 
up with the result 

1 ft nt' /■-|-c» 2 a 

Aiv + Av ^ 2nj-^ dt' di" V / r' dr' ^e'"^^'^">^''-'"^ -^{\G{n,r,r' ,nn')\^)iu{r'). (95) 

''Wo Jo „ Jo *" c'r 

This expression shows that the correlation function appearing under the integral sign only depends on |<' — <"|. Using 
the identity (|5^ we find, for t — >■ +00, that 

1 /•+00 /•+00 o 

Aiv + Av = 2ttj— cZs V / r'dr'n^e"'^^-^'^' — {\G{n,r,r',nn')\'')Lu{r'). (96) 

Using the symmetry s — )■ — s, we can replace /„ ds by (1/2) /__^ ds. Then, using the identities (|65|) and (p6|) . we 
obtain the expression 

A/v + Av = 27r27— y / r'dr'\n\ — (\G(n,r,r',nn)\^)S(n-n')uj(r'). (97) 

r-^ -^ — ' /n or 

Finally, summing Eqs. dM]), dH]) and ([97]), we get 

A// + Am +Aiv+Av = —. (98) 

Therefore, recalling Eq. ([7^. the complete expression of the drift term is 

A^2Tr^-f-y^ f °" dr'\n\\G(n,r,r',nn)\^S{n-n')^ + ^. (99) 

r ^ Jq dr dr 

E. Connection betvi^een the kinetic equation (|45p and the Foliker-Planck equation (|57p 

We have established that the diffusion coefficient and the drift term are given by Eqs. (|55| and ([M|) . Introducing 
the notation (|46|) . they can be written 

1 r^°° 

D = 27727— / r' dr' x(r, r')^(0 - O')w(r'), (100) 

?■ Jo 

and 

A^2Tr^-f- f dr'xir,r')S{n^n')^{r') + ^. (101) 

r Jq dr dr 

Comparing Eq. (|10ip with Eq. ([55| , we see that the "essential" part of the drift term is 

-1 /'+00 J 

V,„ft = 27T'j- dr'x{r,r')S{n-n')^{r'). (102) 



r 



dr 



15 
On the other hand, using an integration by part in the first term of Eq. (|10ip . the total drift can be written 

A = 2.^7/^ rr'dr'u:{r') [-- - -— j x(r,r')^(l^ - O')^- (103) 

Finally, using Eqs. (|100p and (|102p . we find that the Fokkcr-Planck equation ([57)) becomes 

f - 'Ki r^'*'x(--.--'M(o- 0', (vf - y^) , (104, 



When collective effects are neglected, we recover the results obtained in |27h31|| by a different method. As observed 
in our previous works, we note that the form of Eq. ()104p is very similar to the form of Eq. pS]). This shows that the 
Fokker-Planck equation (|104p . with the diffusion coefficient ()100p and the drift term ()102p . can be directly obtained 
from the kinetic equation (|45p by replacing the time dependent distribution a;(r',t) by the static distribution w(r') 
of the bath. This procedure transforms an integro-differential equation (|45|) into a differential equation (|104p |29| . 
Although natural, the rigorous justification of this procedure requires the detailed calculation that we have given here. 
In fact, we can understand this result in the following manner. Equations (|45|) and (|104p describe the evolution of the 
distribution function of a test vortex (described by the coordinate r) interacting with field vortices (described by the 
running coordinate r'). In Eq. (|45p . all the vortices are equivalent so that the distribution of the field vortices u;(r', t) 
changes with time exactly like the distribution of the test vortex a;(r, i). In Eq. (jl04p . the test vortex and the field 
vortices are not equivalent since the field vortices form a "bath" . The field vortices have a steady (given) distribution 
ijj{r') while the distribution of the test vortex w(r, t) = N^P{r, t) changes with time. This distinction was particularly 



visible in our first derivation of the kinetic and Fokker-Planck equations using projection operator methods [28 1 



F. Monotonic profile of angular velocity 

If the profile of angular velocity 51(r) of the field vortices is monotonic, then using the identity 5{Q,{r) — il(r')) 
5{r ~ r')/\D,'{r)\, we find that the expressions of the diffusion coefficient and of the drift simplify into 

2 X{r-,r) 



D{r) = 27r^7^gT-T^(0, (105) 



and 



2 x(r,r) du ^ 



V^-/*=2.^7^^(0, (106) 

where S(r) ~ rQ,'{r) is the local shear created by the field vortices. Comparing Eqs. (|105p and (|106p . we obtain 

V,r^ft = D{r)'^. (107) 

dr 

This relation is valid for an arbitrary distribution of the field vortices provided that the profile of angular momentum 
is monotonic so that their distribution is steady (this important relation was first given in Eq. (123) of Ref. [28J'). 
It can be viewed as a generalization of the Einstein relation (see next section) for a bath that is out-of-equilibrium. 
Therefore, the Fokker-Planck equation ((57| can be written 



dP 1 d 
dt r dr 



fdP dlnoj 



(108) 



This is a drift-diffusion equation describing the evolution of the distribution P{r, t) of the test vortex in an "effective 
potential" Uejf = — Inw(r) produced by the field vortices. From this equation, we find that, for t — > +oo, the 
distribution of the test vortex becomes equal to the distribution of the bath: Peij") = Lu(r)/Nj at equilibrium. 

If we neglect collective effects, we recover the results obtained in |27H3l| . For the usual potential of interaction, the 
function Xbareir^r) is given by (see Appendix [^ : 



1 

8^2 



Xhareir,r) = ^^li-iA, (109) 
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where InA = X]m=i -'^/"^ i^ ^ logarithmically diverging Coulomb factor that has to be regularized appropriately (see 
below). Combining Eqs. (|105p . (|106|) and (|109p . the diffusion coefhcient and the drift term can be rewritten 

„ , , 7 In A , , 7 In A dw , , , , 

^('^)-iM^-W' ^-^^-iM^d^W. (110) 

These expressions were first derived for a thermal bath in [27| and extended to an arbitrary (steady Euler stable) 
distribution of the field vortices in [28|. They also appear in the works |23l426| with a different interpretation (see 
the conclusion for more details). An important feature of these results is that the diffusion coefficient and the 



drift term are inversely proportional to the local shear, a feature first noted in j27| . Furthermore, the diffusion 
coefficient is proportional to the vorticity and the drift term is proportional to the gradient of the vorticity. Since 
"^drift — {D/uj)Vuj in vectorial form, we find that, due to the drift, the test vortex ascends the vorticity gradient. 
This is a purely deterministic effect that appears when the distribution of the field vortices is spatially inhomogeneous 
(in the absence of vorticity gradient, there is no drift). As shown in [23, [2g (see also [2j, [23|), the systematic drift 
of the test vortex is due to a polarization process: the test vortex modifies the distribution of the field vortices and, 
in response, the retroaction of this perturbation causes the drift of the test vortex. In the absence of fiuctuation, 
the test vortex would reach a maximum of background vorticity where Vw = 0. In fact, this systematic effect is 
counterbalanced by the diffusion term that tends to disperse the test vortex. Finally, an equilibrium state results 
from these two antagonistic effects in which the distribution of the test vortex coincides with the distribution of the 
bath. This equilibration process is governed by the Fokker-Planck equation (jlOSp . 

The logarithmic divergence of the Coulomb factor InA (that persists if collective effects are included) was first 
noted in [231 ■ It is due to the failure of the kinetic theory at small scales where collisions between vortices are strong. 
Phen omenologically, the logarithmic divergence can be regularized by introducing cut-offs so that InA = ln(r/d) 
[llll | . A detailed calculation of the lower cut-off d has been made by Dubin & Jin [2^, [2a] . They propose to take 
d ~ Max((5, /) where I is the trapping distance I = (7/47r|E|)^'^ and S is the diffusion-limited maximum separation 
6 = (4£'/|I]|)^/^ where D is the diffusion coefficient given by Eq. (jllOp . Orders of magnitude indicate that r/l ^ 
i?(|E|/7)i/2 - i?(r/7i?2)i/2 _ jsfi/2 (.^vhcrc R is the system size), S/l - (0/^)^/^ - (InA)!/^ ~ [ln(ln7V)]i/2 and 
r/S ^ A^^/^/[ln(lniV)]^/^. Therefore, at leading order, the Coulombian factor scales with N like [38| : 

lnA--ln7V. (Ill) 

2 ^ ' 

This expression may be substituted in Eq. piOp if we are interested in orders of magnitude. 

G. Thermal bath: Boltzmann distribution 

For a thermal bath, the field vortices have the Boltzmann distribution of statistical equilibrium 

Lo{r) = Ae~'^^'^''''''\ (112) 

where ip*{r) = 'ip{r) + (r2L/2)r^ is the relative stream function taking into account the invariance by rotation of the 
system. We note the identity 

^ = -/37c.^ - /?7^f^ - f^L)r, (113) 

dr dr 

where we have used —dip/dr — fir. Substituting this relation in Eq. (|102p . we obtain 

Vd„ft = 2tt^j^(3- / dr'x{r, r')Sin - n'^ifl' - ^lY . (114) 

^ Jo 

Using the (5-function to replace 51' by 17, then using i7 — f2i — —{l/r)d'ip^/dr, and comparing the resulting expression 
with Eq. ^UU\), we finally find that 

Vd„ft = -DP^^. (115) 

In vectorial form, this can be written "V drift = -~D/3^Wip* [27[. The drift is perpendicular to the relative mean field 
velocity u* = — z x Wip* and the drift coefficient (mobility) satisfies an Einstein relation 

C^D/S-f. (116) 
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We note that the drift coefficient ^ and the diffusion coefficient D depend on the position r of the test vortex and 
we recall that the inverse temperature /3 is negative in cases of physical interest. We stress that the Einstein relation 
is valid for the "essential" drift Vdrift, not for the total drift A that has a more complex expression due to the term 
dD/dr. We do not have this subtlety for the usual Brownian motion where the diffusion coefficient is constant. For 
a thermal bath, using Eq. (|115|) . the Fokker-Planck equation ([57]) takes the form 



dP _1 d 
dt r dr 



-«(f-«--t 



(117) 



where D{r) is given by Eq. pOOp with Eq. (|112p . This equation has a form similar to the familiar Smoluchowski 
equation in Brownian theory. This is a drift-diffusion equation describing the evolution of the distribution P(r, t) of 
the test vortex in an "effective potential" C/e// = V'* produced by the field vortices. For t — > +oo, the distribution 
of the test vortex relaxes towards the Boltzmann distribution: Pe{r) = {A/N^)e~^^^''^^^ at equilibrium. Of course, 
if the profile of angular velocity of the Boltzmann distribution is monotonic, we find that Eq. p06p with Eq. (|112p 
returns Eq. (J115I) with a diffusion coefficient given by Eq. (|105l) with Eq. (|112p . However, Eqs. (|115|) and (|100p 
are valid even if the profile of angular velocity in the Boltzmann distribution is non-monotonic. Finally, we note 
that the systematic drift 'Vdrift = —D{r)(3^'V_iju of a point vortex |27| is the counterpart of the dynamical friction 
F/ric = — i?||('y)/3TOV of a star in a cluster [S^, and the Smoluchowski-type Fokker-Planck equation (|117p is the 
counterpart of the Kramers- Chandrasekhar equation (see [3l[ for a development of this analogy). 

H. Gaussian vortex 

For illustration, we consider the particular case of a Gaussian distribution of field vortices 

-M - ^e--V.. (US) 

It can be viewed as a particular statistical equilibrium state of the form pi2p corresponding to /3 — > and ili — > +oo 
in such a way that the product A = P^l remains finite [12l | . It is also similar to the familiar Maxwellian distribution 
of the velocities in a gas at inverse temperature A (provided that r is replaced by v and 7 by m). The Fokker-Planck 
equation (|108p or (|117p describing the relaxation of a test vortex in a Gaussian bath is 



dP I d 
dt r dr 



dP 

rD{r) [ -— + XjPr 



(119) 



where D(r) is given by Eq. (jlOSP with Eq. ()118p . This equation has a form similar to the familiar Kramers equation 
in Brownian theory (provided that r is replaced by v and 7 by m). This is a drift-diffusion equation describing the 
evolution of the distribution P{r,t) of the test vortex in an "effective potential" Ueff = r^/2 (quadratic) produced 
by the field vortices. For t — )• 4-oo, the distribution of the test vortex relaxes towards the Maxwellian distribution of 
the bath: Pe{r) = (A7/27r)e"^'^'''/2 at equilibrium. 

Using fl{r) = r^^ JJ uj{r')r' dr' (see Sec. 6.1 of J3g), the profile of angular velocity of the bath is 

^{r)^:^il-e-^''^'/'). (120) 

If we neglect collective effects, the diffusion coefficient is given by Eq. (|110p . Using Eqs. (jllSp and (|120p . we obtain 
the expression 

/ N 7^AlnAr^ 

^('^) - 4(2eJv.-A7.^-2) - ^'''^ 

We first note that the diffusion coefficient diverges like D{r) ^ In A/(Ar^) for r — > due to the vanishing of the shear 
in the core of the vortex: I](r) ~ — (A^7^r/87r)r^ as r — )• 0. This indicates a failure of the kinetic theory for r — >■ 0. 
In fact, the kinetic theory developed in this paper is valid for sufficiently large shears. In the absence of shear, the 
expression of the diffusion coefficient is different (and finite) as discussed in [42,173]. Therefore, the expression (|12ip 
of the diffusion coefficient is valid for sufficiently large r. We note that the diffusion coefficient decays very rapidly at 
large distances since D{r) ^ {'j'^X\nA/8)r^e~^''^ 1"^ for r -^ +00. More generally, using Eq. (J105I) and the results of 
Sec. 6.1. of [38|, we have P>(r) ^ (27r'^7/r)x(?', ?')w(r)r^ for r — > +cx) for an arbitrary distribution a;(r) of the field 
vortices. 
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The Fokkcr-Planck equation (|119[) with the diffusion coefficient (|12ip has been studied in [S^ by applying results 
previously obtained in the context of tlie HMF model. For this model, the equivalent of the Fokker-Planck equation 
(|119p with Eq. (|121|) . where the position r is replaced by the velocity v, has been studied in [4a,|4^. This Fokker- 
Planck equation presents unusual features because the diffusion coefficient decreases very rapidly with the distance. 
This leads to "anomalies" with respect to the usual Brownian motion. By applying the approach of Bouchet & 
Dauxois [46|, one finds [381 that the auto-correlation function of the test vortex (r(O)r(r)) decays algebraically like 
Ini/i (this algebraic decay was first obtained by Marksteiner et al. [7J1 for the logarithmic Fokker-Planck equation 
to which Eq. (J119p can be mapped). On the other hand, by applying the approach of Chavanis & Lemou |48|, one 
finds [S^ that the normalized distribution u{r, t) — P{r, t)/Pe{r) has a front structure and that the front evolves very 
slowly with time, scaling like rf{t) oc (Ini)^/^ for t — > -|-oo. These results show that the relaxation of Pjr, tVtowards 
the equilibrium distribution Pe{r) is not exponential. This is intrinsically due to the absence of gap [3^|4y, [tJ in the 
spectrum of the Fokker-Planck equation (|119p with diffusion coefficient (|12ip . 

I. The relaxation time of a test vortex in a bath 

The above derivation of the Fokker-Planck equation, relying on a bath approximation, assumes that the distribution 
of the field vortices is "frozen" so that their vorticity profile Lu(r) does not evolve in time. This is always true for a 
thermal bath (J112p . corresponding to a distribution at statistical equilibrium (Boltzmann), because it does not evolve 
at all. For a stable steady state of the 2D Euler equation with a monotonic profile of angular velocity, this is true 
only on a timescale shorter than the relaxation time tfj of the system as a whole. However, as we have indicated in 
Section IIIIBi this timescale is very long because the relaxation time tR of the system as a whole is larger than Ntn ■ 

Recalling that 7 - 1/iV, the Fokker-Planck operator in Eqs. !(TUE\\ and PTT)) scales like \nN/N (see Eq. pTTj) for 
the logarithmic correction). Therefore, the distribution P{r, t) of the test vortex relaxes towards the distribution uj{r) 
of the bath on a typical time 

tT" - ^tD. (122) 

This is the timescale controlling the relaxation of the test vortex, i.e. the time needed by the test vortex to acquire the 
distribution of the bath (see Appendix [Dl for a more precise estimate). It should not be confused with the timescale 
(j48| controlling the relaxation of the system as a whole. Since the timescale i^*'' is shorter than the timescale Ir 
on which uj{r,t) changes due to collisions (in the case where uj{r,t) has a monotonic profile of angular velocity), 
we can consider that the distribution of the field vortices is "frozen" on the timescale (|122p . Therefore, our bath 
approximation is justified on this timescale and our approach is self-consistent. 

V. CONCLUSION 

In this paper, we have developed the kinetic theory of point vortices by taking collective effects into account. We 
have used the formalism of Dubin h O'Neil [20| valid for axisymmetric mean flows. This improves our previous works 
[27h31| where these effects were neglected. Collective effects amount to replacing the "bare" potential of interaction 
by a "dressed" potential, without alte ring the overall structure of the kinetic equation. In plasma physics, collective 
effects included in the Lenard-Balescu [2l|, l22| equation are important because they take into account screening effects 
and regularize, at the scale of the Debye length, the logarithmic divergence at large scales that appears in the Landau 



37[ equation. In the case of point vortices, there is no divergence in the Landau-type kinetic equation (|43p that ignores 
collective effects, so that their influence may be less important than in plasma physics. We have also developed a 
test vortex approach and a Fokker-Planck theory. We have obtained the expressions of the diffusion coefficient and 
drift term directly from the equations of motion, taking collective effects into account. We have shown that they can 
also be obtained from the kinetic equation by making a bath approximation leading to a Fokker-Planck equation. We 
have presented the results for axisymmetric flows, but similar results can be obtained for unidirectional flows |27l.l28l|. 
The kinetic theory developed in this paper is valid at the order 1/iV so that it describes the evolution of the system 
on a timescale Ntu. This is sufficient to study the relaxation of a test vortex in a bath of field vortices since the 
corresponding relaxation time is of order {N/ lnN)t£) (see Section HVip . However, this is not sufficient to describe the 
evolution of the system as a whole towards the statistical equilibrium state (Boltzmann) because, for axisymmetric 
flows, the relaxation time is larger than Ntjy (see Section UlIBp . Therefore, we need to develop the kinetic theory at 
higher orders (taking into account three-body, four-body,... correlations). At present, this is not done and the scaling 
of the relaxation time with A^ remains an open problem. 



19 

We can try to make speculations by using analogies with other systems with long-range interactions which present 
similar features [3l|. The most natural scaling would be N'^tjj which corresponds to the next order term in the 
expansion of the basic equations of the kinetic theory in powers of 1/N [2^, |3l|, |6J| . An N^ scaling is indeed obtained 
numerically |75l . |76| for spatially homogeneous one-dimensional plasmas for which the Lenard-Balescu collision term 
vanishes at the order 1/N |5ll. |52||. However, the scaling of the relaxation time may be more complex. For example, 
for the permanently spatially homogeneous HMF model (for which the Lenard-Balescu collision term also vanishes 
at the order 1/A^ (ia. |47| ) . Campa et al. [73| report a relaxation time scaling like e'^to (this timescale is, however, 
questioned in recent works [78| who find a N'^to scaling). It could also happen that the point vortex gas never 
reaches statistical equilibrium, i.e. the evolution may be non-ergodic [l^. All these speculations [31| could be checked 
numerically by solving the iV- vortex dynamics. The situation should be different for more general flows that are 
not axisymmetric. In that case, there are potentially more resonances so that the relaxation time could be reduced 
and achieve the natural scaling Ntjj corresponding to the first order term in the kinetic theory [2^, [31| . Similarly, 
one dimensional spatially homogeneous systems with long-raiige interaction are expected to relax towards statistical 
equilibrium on a timescale Nto due to additional resonances |53ll . An Nto scaling is indeed obtained numerically for 
spatially inhomogeneous one dimensional stellar systems J79l - l82l |. for the spatially inhomogeneous HMF model [83] and 
for the relaxation of a non-axisymmetric distribution of point vortices |7l| . On the other hand, for the HMF model, 
if an initially spatially homogeneous distribution function becomes Vlasov unstable during the collisional evolution, a 
dynamical phase transition from a non- magnetized to a magnetized state takes place (as theoretically studied in |84{ ) 
and the relaxation time could be intermediate between N'^t^ (permanently homogeneous) and Ntjj (permanently 
inhomogeneous). In that situation, Yamaguchi et al. [83] find a relaxation time scaling like N^tn with S = 1.7. 
The previous argument (leading to 1 < 5 < 2) may provide a first step towards the explanation of this anomalous 
exponent. The same phenomenon (loss of Euler stability due to collisions and dynamical phase transition from an 
axisymmetric distribution to a non-axisymmetric distribution) could happen for the point vortex system. 

Finally, we would like to conclude by briefly mentioning other works on the kinetic theory of point vortices that 
are related to our study. Schecter & Dubin [2j, [2J] have considered the motion of a point vortex in a background 
vorticity gradient (without fluctuation) . They obtain an expression of the drift that coincides with the one obtained 
previously in [27ll28{. However, the point of view is different because they consider the evolution of a test vortex in an 
external vorticity field without fluctuation, while we consider the evolution of a test vortex in a bath of field vortices 
that creates a vorticity gradient but also contains fluctuations. The vorticity gradient leads to a (deterministic) drift 
as a result of a polarization process and the fluctuations lead to a diffusion. Therefore, the test vortex has a stochastic 
motion that can be modeled by a Fokker-Planck equation of a drift-diffusion type. There is no Fokker-Planck equation 
in the approach of Schecter & Dubin [2J, [2J] since the evolution of the test vortex is purely deterministic. On the other 
hand, Dubin & Jin |25| consider the diffusion of point vortices in an external shear. They obtain an expression of the 
diffusion coefficient that coincides with the one obtained previously in [23, UM from the Kubo formula (in particular, 
they recover the shear reduction found in [27[). However, the point of view is different because they consider the 
effect of an external shear without vorticity gradient. Therefore, there is no drift. By contrast, in our approach, the 
test vortex moves in a bath of field vortices that presents a vorticity gradient. In addition to its diffusive motion, 
the test vortex modifies the distribution of the field vortices (as in a polarization process) and the retroaction of the 
field vortices causes the drift of the test vortex. In a sense, the works of Dubin and collaborators isolate the effects 



of drift |23l |2J] and diffusion [25| while they arc taken into account simultaneously in our approaches |27l - [31| . The 
extension of the kinetic theory to the case of point vortices with different values of individual circulation 7 has been 
performed in [2^, [S^, ISGJ - ISSf . The statistics of the velocity fluctuation arising from a random distribution of point 
vortices has been studied in [TJ, I89l - l9ll | by analogy with the study of Chandrasekhar & von Neumann [O^] on the 
statistics of the gravitational force arising from a random distribution of stars. This theory has been used in [73| to 
obtain an expression D ex 7(ln A^)^/^ of the diffusion coefficient of point vortices in the absence of shear (i.e. in the 
opposite limit to the one considered here). This result has been applied to the problem of 2D decaying turbulence in 
[93| in order to obtain the expression of the exponent of anomalous diffusion (x^) ^ t^+?/^ where ^ is the exponent 
controlling the decay of the vortex number A^ ^ t~^ (it is argued in [9j, [9J] that ^ = 1 in the strict asymptotic scaling 
regime t — > -l-oo so that v = 1 + ^/2 = 3/2). The energy spectrum of the point vortex gas has been determined in 
(43. lOll. I95I - I97I I . Finally, a quasilinear theory of the 2D Euler equation has been developed in |98| in orde r to describe 



a phase of "gentle" collisionless relaxation. For general references about vortex dynamics, see [10d . ll01| . 

Appendix A: Potential of interaction 

We assume that the stream function ip{r) is related to the vorticity Lu{r) by an equation of the form Lip = —u, 
where L is a linear operator of the Laplacian type. The potential of interaction u(r,r'), which is the Green function 
of the operator L, is defined by Lw(r,r') = —S{r — r'). For simplicity, we shall consider an infinite domain where 
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u = u{\r — r'\) only depends on the absolute distance between two points (our results can be extended to more general 
situations by using the Green functions of Lin [93 )• In that case, the stream function is related to the vorticity field 
by an expression of the form 

ij{r,t)= u{\r-r'\)uj{r',t)dr. (Al) 

It can be written tjj = u * uj where * denotes the product of convolution. For ordinary flows, the stream function is 
related to the vorticity by the Poisson equation /S.ip = — w. In that case i = A and the potential of interaction in 
an infinite domain is given by u(|r — r'|) = — (l/27r)ln|r — r'| which corresponds to a Newtonian (or Coulombian) 
interaction in two dimensions. In the quasigeostrophic (QG) model describing geophysical flows 140|, the stream 



function is related to the (potential) vorticity by the screened Poisson equation A-0 — kj^tp 



where kr, is 



the Rossby radius. In that case L = A — fc|j and the potential of interaction in an infinite domain is given by 
u(\y — r'l) = (l/27r)A'o(fcfl,|r — r'|). The previous results are recovered for fc/j — s- 0. 

Let us first consider the ordinary potential corresponding to L = A. Introducing polar coordinates, the Poisson 
equation can be written 



I d d5i} 1 d'^Sip 
r dr dr r^ 89^ 

Taking the Fourier-Laplace transform of this equation, we obtain 



-Sio. 



1 d d n? 
r dr dr r"^ 



Sip{n, r, <t) = —5Lb{n^ r, cr). 



(A2) 



(A3) 



This is of the form C5^ 



-(5a) with £= -|-r|- 

r or or 



Substituting Eq. ^ in Eq. (|M1), we find that 



1 d d 



r or 

r dr dr r^ nfl — a 



S'ip{n,r,a) 



6ili{n, r, 0) 
i{nn — a) 



(A4) 



Therefore Stp{n,r,a) is related to (5(i(n, r, 0) by a relation of the form (P5)) where the Green function G{n,r,r' ,a) is 
defined by 



1 d d 

r dr dr 



nfl 



G{n,r,r',a) 



S{r ~ r') 
lixr 



If we neglect collective effects in the foregoing equation, we obtain 



i_a d_ 

r dr dr 



Gbare{n,r,r') 



5{r - r') 
2'nr 



(A5) 



(A6) 



This shows that the bare Green function Gtarei'iT'i ^i f') — Uni^i ^') is the Fourier transform of the potential of inter- 
action u that is solution of the Poisson equation Au = —5. 

Let us now consider a general form of interaction between point vortices. The fluctuations of the stream function 
are related to the fluctuations of the vorticity by 



5i}{r,t)= I u{\T-v'\)5uj{v',t)dT. 
The potential of interaction can be written 

u(|r - r'l) = u (yjr"^ + r'^ - 2rr' cos{6 - 9')) = u{r, r', 0), 
where cj) — 9 — 9' . Due to its ^-periodicity, it can be decomposed in Fourier series of the form 

/>27r II 

M(r, r',0) = y^e™'^M„(r, r'), u„(r, r')= / --u{r,r' ,<f)) cosincj)). 



(A7) 



(A8) 



(A9) 
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Taking the Fourier-Laplace transform of Eq. (|A7p and using the fact that the integral is a product of convolution, we 
get 



r+oo 

Sip^n, r, a) = 2Tr / r'dr' ■u,n{r,r')6u]{n,r' ,a). 
Jo 

If we introduce Eq. ([23|) in Eq. (JA10| , we obtain a Fredholm integral equation 

S^Pin,r,a)+2n r'dr' u„ir,r')^^S^in,r' , a) ^ 27r r'dr' Mr,r') .^ ^^] ' ' 

Jo ml' — a Jq i(ml' ~ a) 

which is equivalent to Eq. ([M|) . If we neglect collective effects, the foregoing equation reduces to 



5il){n,r,a) =271 I r'dr' Un{r,r') ., ' ' ■ 

i[ml — a) 



(AlO) 



(All) 



(A12) 



Comparing Eq. (|A12|) with Eq. (|25p . we see that Gbare{n, r, r') = Un{r, r') is just the Fourier transform of the "bare" 
potential of interaction u. This is equivalent to CGbare{'n,r,r') = —S{r — r')/2Tir. When collective effects are taken 
into account, Eq. (jA12[) is replaced by Eq. ([25)1 where the "dressed" potential of interaction satisfies 



G{n,r,r',(T) +2tt / r"dr" Un{r, r") !"" ^''" G{n, r" , r' , cr) = Un{r, r'), 
Jo nil" - a 



(A13) 



which is equivalent to Eq. ([25]) . 

For the Newtonian interaction u(|r — r'|) = — (l/27r) ln|r — r'|, we have u{r,r' , (j>) = — (l/47r) ln(r^ + r'^ — 2rr'cos( 
The integrals in Eq. (|A9I) can be performed analytically 3lJ, and we obtain 



Therefore 



Un{r,r') 



u(r,r', 



1 fr^ 
47r|n| \7'> 



l"l 



uoir,r') = ---lnr>. 
Ztt 






Jn4> 



(A14) 



(A15) 



which is the Fourier decomposition of the logarithm in two dimensions. The function defined by Eq. (|46p can be 
written [H: 



^ +CXD ^ / \ 2n w 

n— 1 ^ ^ 



l-( — 



Taking r' = r in the foregoing expression, we obtain 

1 ^°° 1 1 

Xbare{r,r) = —Y,- = —\nA, 



n— 1 



(A16) 



(A17) 



where In A is the Coulomb factor (see Sec. lIVFp . Analogous expressions, valid in a circular domain of size R, are 
given in j29| . 



Appendix B: Auto-correlation of the fluctuations of the one- particle density 

According to Eq. (|22p . we have 



27r 



de 



'^"^ de' 



(<5c2;(n,r,0)te(n',r',0)) = ^ ^e~<-'+"'''HSLoi9,r,O)Su:{0' ,r' ,0)) 



27r 7o 27r 



^ / ^e-("«+"'«')[(.;,(0,r,O)c.,(0',/,O))-^(rV(r'))]. (Bl) 



27r 7o 27r 
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The expression (|4]) of the discrete vorticity distribution leads to 



+ ,=5;(,(9-,.)«i_iJ,(9.-<,,):!(rl;iii)) 



= 2TT-foj{r)S{9-e')^^^^—^+oj{r)uj{r'), (B2) 

ZTrr 

where we have assumed that there is no correlation initially (if there are initial correlations, it can be shown that 
they are washed out rapidly [a, [5a] so that they have no effect on the final form of the collision term). Combining 
Eqs. dHl]) and US]), we obtain 

{5uj{n,r,0)6uj{n' ,r' ,0)) ^ -fSn.-n'%-^uj{r). (B3) 

Iirr 

Appendix C: The function K{r,r',t) 

The kinetic equation (j45p can be rewritten in the more compact form 

duj 1 d f+°° , . ^, / ^ A 9 I d \ , , , , , 



(CI) 



with 

K{r, r', t) = 27r^^^x{r, r\t)5{n{r, t) - n{r' , t)). (C2) 

In the thermal bath approximation, the diffusion coefficient and the drift terms can be rewritten 

'•+00 ^^' 



1 r°° r°° dr' duj 

D{r) = - dr' Kir,r'Mr'), Vd„ft{r) = ^ K{r,r')^{r'), 



(C3) 



Air)-rrr;r'.ir'){l^-L^yS^. (C4) 

Jo \r or r or J rr 

Appendix D: A more precise estimate of the relaixation time 

When the field vortices have the Gaussian distribution (|118p . a more precise estimate of the relaxation time can be 
given. We first define the "vortex size" R hy R^ = (r'^) = L/T where L is the angular momentum. For a Gaussian 
distribution, R^ = 2/A7. If we set x = yX'-f/2r, the Fokker-Planck equation (|119l) can be rewritten 

(Dl) 



dP _ 1 1 a 

dt tu X dx 
with 



dP 

xGix) \^— + 2Px 
ox 



^^^) = ;^^-3^3T' ^«=vi^- (D2) 

The "reference time" Ib, gives an estimate of the relaxation time of the test vortex in a Gaussian bath. It we define 
the dynamical time by to = R^ /T = 2/\^T, we obtain 

4r 8iV , , 
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where we have used Eq. (|llip to get the equivalent for N -^ +00. We can also estimate the relaxation time as follows. 
The typical position of the test vortex increases like (r^) ~ AD{R)t. The relaxation time t^ is the time needed by the 
vortex to diffuse over a distance R. Taking (r^) = R^ in the foregoing formula, we obtain 

where we have used G{1) = l/(e — 2) ~ 1.39221. Finally, although the relaxation towards the Gaussian distribution is 
not exponential, a measure of the "relaxation time" is provided by i', = ^^^ where ^ = D{R)Xj is the drift coefficient 
(mobility) given by an Einstein formula. This yields 

1 4 

^'- " DiR)X-f = A72G(l)lnA = ^^"^ ^°^^ 

The same arguments can be extended to the Fokker-Planck equations ()108p and ()117p , leading to the scaling (|122p of 
the relaxation time, up to a numerical factor. 
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